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Abstract 



We provide an extension of the method of asymptotic decompositions of vector fields with finite- 
time singularities by applying the central extension technique of Poincare to the dominant part 

> ■ 

00 ■ of the vector field on approach to the singularity. This leads to a bundle of fan-out asymptotic 

original system. We show how this method can be useful to describe a single-fluid isotropic 
universe at the time of maximum expansion, and discuss possible relations of our results to 
£f~) • structural stability and non-compact phase spaces. 
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1 Introduction 

A popular method to represent infinity in general relativity is Penrose's conformal method (cf. Ref. pQ) 
wherein the overall structure of a physical spacetime is conformally changed so that its infinitely remote 
regions become the boundary of a new, unphysical spacetime. Infinity is then classified according to the 
behaviour of the various geodesies of the new spacetime near its boundary. Until recently, the conformal 
method was the only useful method to study the 'structure of infinity' in relativistic field theories. In finite 
dimensions, the corresponding idea to analyze the global structure of solutions especially for planar systems 
stereographically is due to Bcndixson, cf. Ref. [2]. 
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The method of asymptotic splittings (cf. [3] and refs. therein) offers another asymptotic representation 
of solutions to the dynamical equations near infinity (that is where some component of the solution diverges), 
and it is possible to describe it in three main steps as follows. First, we recall some simple definitions. We 
shall work on open subsets or R™, although our results can be extended without difficulty to any differentiable 
manifold MP. We shall use interchangeably the terms vector field / : M n — > TM n and dynamical system 
defined by / on M n , x = f(x), with (•) = d/dt. Also, we will use the terms 'integral curve' x(t,xo) of the 
vector field / with initial condition x(0) = xq, and 'solution' of the associated dynamical system x(t;xo) 
passing through the point xq, with identical meanings. We say that the system x = f(x) (equivalently, the 
vector field /) has a finite-time singularity if there exists a t s G K and a Xq 6 M n such that for all M G R 
there exists an S > such that \\x(t; Xo)\\lp > M, for \t — t s \ < 6. Here x : (0, b) — > M n , x = x(to) for some 
to £ (0,6), and || ■ \ \lp is any LP norm. We say that the vector field has a future (resp. past) singularity if 
t s > to (resp. t s < to). Note also, that to is an arbitrary point in the domain (0,6) and may be taken to 
mean 'now'. Alternatively, we may set r = t — t s , r G (— t s , b — t*), and consider the solution in terms of the 
new time variable r, x(t; xq), with a finite-time singularity at r = 0. We see that for a vector field to have 
a finite-time singularity there must be at least one integral curve passing through the point xq of M n such 
that at least one of its LP norms diverges at t = t s , and we write 

fim ||x(t;xo)||Lp = oo, (1.1) 

to denote the behaviour of the solution on approach to the finite-time singularity at t s . 

One of the most interesting problems in the theory of singularities of vector fields is to find the structure 
of the set of points xq in MP 1 such that, when evolved through the dynamical system defined by the vector 
field, the integral curve of / passing through a point in that set satisfies property ()1.1|) . The evolutionary, 
geometric nature of this set especially near singularities of the field, is of prime interest in this paper. In 
the method of asymptotic splittings, we imagine that on approach to a finite time singularity, / decomposes 
into a dominant part / and another, subdominant part (consisting possibly of more than one additive 
components f^\ix = 1, • • ■ ,k), as follows: 

/=/ (0) +/ (l) + ... + / « (L2) 

This asymptotic decomposition is highly non-unique. Then, in the first step, for each particular decomposi- 
tion we drop all terms after /W and replace the exact equation x = f(x) by the asymptotic equation 

x = fM(x). (1.3) 

Repeating this for all possible decompositions of /, we end up with a bundle of asymptotic dynamical 
systems to work with. In Step two, we look for scale-invariant solutions (called dominant balances) in each 
asymptotic system. Each such system may contain many possible dominant balances resulting in various 
possible dominant behaviours near the singularity. Balancing each asymptotic system, requires a careful 
asymptotic analysis of the subdominant parts in the various vector field decompositions. In the last step, we 
check the overall consistency of the approximation scheme and build asymptotic solutions for each acceptable 
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balance in a term-by-term iteration procedure that ends up with a formal series expansion representation of 
the solutions. When the whole procedure is completed, we have finally constructed formal series developments 
of particular or general solutions of the original system of equations, valid in a local neighborhood of the 
finite-time singularity. From such expansions we can deduce all possible dominant modes of approach of the 
field to the singularity, decide on the generality of the constructed solutions and also determine the size and 
part of the space of initial data that led to such a solution. However, this method has two shortcomings: 

1. It does not give any information about the qualitative behaviour of individual orbits 

2. It does not distinguish between the behaviour of those solution components tending to +00 from those 
ones tending to — 00. 

The qualitative method of central projections and compactifications due to Poincare is usually presented, 
cf. Refs. [H O [6l [7] , as suitable only for homogeneous vector fields and singularities at infinity, that is when 
we have a homogeneous vector field diverging as t — > 00. The purpose of this paper is to emphasize the fact 
that this method, when considered from the viewpoint of asymptotic splittings, is really suitable for any 
weight-homogeneous vector field that blows up at a finite time (or at infinity) . In this way, we may suitably 
adapt the ingenious method of Poicare's to study the finite-time singularities that arise in cosmology and 
relativistic field theory, wherein a healthy dose of problems having precisely this kind of vector field and 
asymptotic structure appear with phase spaces of any finite number of dimensions. It is hoped that this will 
provide a useful alternative to the stereographic techniques of Bendixson and Penrose mentioned above. 

In this paper, we revisit the method of central projections and use it as a suitable complement and 
development of the method of asymptotic splittings, in order to decide about the qualitative behaviours of 
the asymptotic orbits that result from the various possible asymptotic decompositions. In the next section, 
we give the details of the method of central projection for one-dimensional systems, and describe the state 
of maximum expansion for a universe with a single fluid in general relativity, a problem known not to admit 
a complete description with standard methods. In Section 3, we describe the amalgamation of asymptotic 
splittings and central projections in two dimensions, and in the last section we discuss various issues that lie 
ahead and result from the approach taken in this paper. 

2 The one-dimensional case 

Let us consider the one-dimensional system x — f(x) and assume that there is a finite time t s such that the 
solution x has a singularity, x(t s ) = +00, or x(t s ) = —00. In general relativity (which is the principal area of 
application we have in mind for the techniques discussed here), one is accustomed to view such singularities 
as following from other plausible physical and geometric hypotheses - the so-called singularity theorems, cf. 
[9j [lOl QT]. We imagine the one-dimensional phase space T = {x : x G R} of the system as sitting at the 
(X, 1) (cotangent) line of an (X, Z) plane, and we let x be any particular state of the system on this line (so 
that x = X in this case). The Poincare central projection pr c is a map from T to the upper semi-circle S\_ 
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sending each phase point to a point on the circle such that 

pr c : T -> Si : x i-> 0, with a; = cot6». (2.1) 

In this way, the two possible infinities of x at the two ends of T are bijectively separated (note that this 
cannot happen with a stereographic map like Penrose's), 

x = +oo -)■ 6 = (2.2) 
x = -oo -> = 7T, (2.3) 

the new phase space, the so-called Poincare's circle, <Si is now compact, # S [0, n], and the original system 
x = f(x) reads, 

9 = -sin 2 (9/(cot 9) = g(9). (2.4) 

So past singularities are met as 9 \, 0, while future ones are at 9 f 7r. Let us suppose for the sake of 
illustration that we are interested in the asymptotic properties of the solution near a past singularity. In 
the next step of the Poincare central projection method, we are after a singular asymptotic system on the 
compactified phase space 5+. According to the method of asymptotic splittings, the field / is asymptotically 
decomposed into two parts, a dominant one /(°), and another subdominant, /( sub \ i.e., / = /(°) + /( sub ), 
in a highly non-unique manner. We take any particular decomposition for which = ax M ,M S Q, and 
asymptotically as 9 \. we find that g{9) = —a9 2 ~ M , that is we arrive at the singular asymptotic system 

6 = -a9 2 - M . (2.5) 

As in the method of asymptotic splittings, we need to repeat the whole procedure below for all possible 
decompositions. Notice that this system is only valid locally around the singularity at 9 = 0, not everywhere 
on the compactified phase space 5i. The last step is now to obtain a singularity-free system valid locally 
around the singularity. For this purpose, one changes the time t to a new one, r, given by d/dt — 9 X ~ M d/dr, 
to arrive at the complete asymptotic system defined on the Poincare circle, 

9 = -a9, (2.6) 

and valid near 9 = 0. We conclude that the system has a past attractor when a > 0, otherwise all orbits are 
repelled near the singularity. 

To assess the merits of this method, we revisit the evolution of the single-fluid FL models in general 
relativity with a linear equation of state p = (7 — l)p, studied in Ref. [5], pp. 58-60. Introducing the density 
parameter Q defined in terms of the fluid density p and Hubble rate H by ft = p/3H 2 , the evolution of the 
system is governed by the equation 

n = -(37-2)(i-n)n, (2.7) 

and it is well known that for closed models, after a finite time interval, at the instant of maximum expansion, 
fl blows up to +00, thus making the description of the whole evolution incomplete (cf. Ref. jS], p. 60). To 
describe this, we choose the decomposition having f(°> = ail 2 , where a = 3j — 2 and M = 2. The resulting 
asymptotic system is 9 = —(37 — 2), and taking 9 = at t = t m3lX , we find that 9 = —(37 — 2)(t — i m ax)- 
The central projection gives f2 = cot 9 = cot(— (37 — 2)(t — t max )), near the point of maximum expansion. 
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3 Dimension 2 



For homogeneous vector fields, the various reductions involved in central projections are described in Refs. 
[IJ EJ [7] and will not be repeated here. Below, as in the method of asymptotic splittings, we focus 
on weight- homogeneous vector fields and finite-time singularities exclusively. Any C 1 function h is called 
weight-homogeneous if we can find a d G K and a nonzero vector w € W 1 such that for all numbers t, we have 

h(t Wl Xl ,t W2 x 2 ,--- ,t w "x n )=t d h{ Xl ,x 2 ,--- ,x n ). (3.1) 

We then speak of the weight vector w — (wi) and the weighted degree deg(h, w) = d of the function h. For 
instance, the function fx = 2xy is weight-homogeneous with weight vector w = (—1,-1) and weighted degree 
d = —2, and /a = 4y 2 — x 2 is also weight-homogeneous with the same w and d. When w = (1, • • • ,1), we 
say that ft. is a homogeneous function. 

A weight-homogeneous vector field f = (fx, ■ ■ ■ , f n ) is one for which we can find a a common weight vector 
w for the f^s. When the components = 1, • • • , n, have weighted degrees di = deg(/i, w),i = 1, ■ • • , n„ 
we say that the vector field / has weighted degree deg(/, w) — d with d — (di). For example, the vector field 
/ = (4y — x 2 /2, xy + x 3 ) has weight vector w = (—1, —2) and weighted degree d — (—2, —3). 

An important subclass of the weight-homogeneous vector fields are the scale-invariant vector fields, for 
which their weight vectors and degrees are related by, 

d = w-l, (3.2) 

with 1 here meaning the vector having 1 in every slot (hence, = Wi — 1). Thus the degree d of a scale- 
invariant vector field measures its failure from being exactly homogeneous. An obvious property of any 
scale- invariant field / with weight vector p, is that the 'scale-invariant' system x = f admits a solution of the 
form x = at p , a £ C™ (in this notation, Xi — ait Vi in components), when the algebraic system a-p = f(a) has 
a nontrivial solution (that is some component of a is nonzero). In this case, the solution x = at p is invariant 
under the scaling t — > Xt,x — > X p x. 

In the method of asymptotic splittings, we look for weight-homogeneous decompositions, that is vector 
field decompositions of the form (|1.2[) such that the following conditions hold: 

1. all vector field 'components' / W in the sum are weight- homogeneous with the same weight vector w, 

2. their weighted degrees are given by d^' = w — j — 1, and d^ > d^ > ■ ■ ■ > d^ k \ 

3. the most nonlinear, additive component /^°^ is scale-invariant. 

For the sake of illustration, we shall confine our attention to planar systems of the form 

x = P(x,y), y = Q(x,y), (3.3) 

and we shall assume that this system has a finite-time singularity at time i s , as in the definition (jl.ip 
(everything we do generalizes easily to any finite dimension). As in the 1-dimensional case of the previous 
section, since the original system (|3.3[) is incomplete, we are eventually after a complete, 'asymptotic' system 
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defined on the 2-dimensional sphere S 2 . Taking the phanar phase space (x,y) to be the Z = 1 plane of a 
new 3-space with coordinates X, Y, Z, the embedding 

pr c : R 2 ->• Si : (x, y) ^ {X, Y, Z) : x = X/Z, y = Y/Z, X 2 + Y 2 + Z 2 = l, (3.4) 

interprets any blow up in the solution vector (x,y), x,y — >• oo, as the limit Z — > 0. (We note that any sign 
in the x,y's is taken care of by the signs of the new X, Y variables, and so we can limit ourselves on the 
positive semi-sphere, that is consider only the direction Z J. 0.) Thus infinity in x, y is now (at) the Z = 
circle of infinity, X 2 + Y 2 = 1. This is the first step of the so-called Poincare compactification for the system 
(|3.3[) wherein the embedding pr c is the central projection transform, and sends a planar phase point (x, y) to 
a point of the semi-sphere (the 'Poincare sphere') S\. Under the map (|3.4p . the system (|3.3[) becomes the 
singular system [U [6j [7] , 

X = Z((l - X 2 )P - XYQ) 

Y = Z(-XYP + (1 - Y 2 )Q) (3.5) 
Z = -Z 2 (XP + YQ). 

We next aim to obtain asymptotically, through a series of transforms, first a reduction of the system (|3.5[) 
that will be complete and valid on 57., and secondly another reduction valid precisely only along the circle 
of infinity. We suppose we have used the method of asymptotic splittings and ended up with the spectrum 
of all possible asymptotic decompositions of the form (|1.2j) of the vector field (P(x,y),Q(x,y)). We shall 
work with a particular weight-homogeneous decomposition that has dominant part given by 

/W =(P(°)Q(°)), (3.6) 

and we suppose that a = (ai),i = 1,2, is the weighted degree of f(°\ We then define, 

M = max{|ai|,i= 1,2}, (3.7) 

and using this, we can introduce the time transform t — > r{t) (monotone when Z > 0), 

so that the singular system (|3.5p becomes the following system for the dominant part of the vector field 
(P,Q) ((') = d/dr), 

Intermediate system: 

x' = z M ((i-i 2 )p(°)-xyQ(°)) 

y' = z M (-xrp (0) + (i-r 2 )Q (0) ) (3.9) 

Z' = -Z M+1 [XP^ +YQ<®y 

There is an important difference with respect to the homogeneous case treated in Refs. [HO 0(7], in that we 
consider only the dominant part of the vector field asymptotically, and this will be different for each particular 



asymptotic decomposition of the field. The 'subdominant' part is truly subdominant asymptotically because 
of the validity of the subdominance condition, cf. [3] , 

lim /S " b( f P) = 0. (3.10) 

In a sense, the standard presentation of the Poincare method corresponds to the all-terms-dominant decom- 
position presently (the only possible one when the field is homogeneous). We can further make this system a 
complete asymptotic system valid on S\ by introducing the sphere vector field {Pg 2 , Qs£ Y, Z) defined 
by the forms, 



P sl (X,Y,Z) = Z"PW(|,|) 



(3-11) 



Q sl (X,Y,Z) = Z M Q(°>(|,|). (3.12) 

Then the intermediate system above becomes one valid on <S? without singularities, namely, 
Complete asymptotic system on S+: 

X' = (1 - X 2 )P sl - XYQ S 



Y' = -XYP sl + (l-Y 2 )Q si (3.13) 
Z' = -Z(XP SI +YQ SI 



We note that in this system, any additive terms proportional to Z are asymptotically negligible on approach 
to the singularity, and so this system asymptotically becomes one valid on the (X, K)-plane along the circle 
of infinity, S 1 : X 2 + Y 2 = 1, that is we have 

Asymptotic system on circle of infinity: 

X 1 = -Y(XQ s i -YP s i) (3.14) 
Y' = X{XQ s x -YP s x), on the circle X 2 + Y 2 = l. (3.15) 

Here we have introduced the circle vector field (P51 , Q51 ) (X, Y) given by 

P s i{X,Y) = P 5 2(X,y,0) (3.16) 
Q s i{X,Y) = Q sl (X,Y,0). (3.17) 

Some remarks about the system (|3.14H3.15"|) defined on the Poincare circle of infinity are in order. First, 
through the series of transforms defined above, for any weight-homogeneous vector field / with a finite-time 
singularity we have found that infinity in the original variables x, y is now an invariant manifold (circle), 
with the dynamics near infinity determined by the constrained 2-dimensional system (|3. 14113. ll)]) . effectively 
1-dimensional, on the equator. Secondly, for any given system of the form p.3[) . there correspond precisely 
the same number of systems of the form (13. 14113.15]) as the number of admissible asymptotic decompositions 
of the original vector field p.3j) . We note that an asymptotic splitting is admissible provided we come up 
eventually with a valid asymptotic solution in the form of a formal series. In particular, we have to find the 
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Kowalevskaya exponents of the various dominant balances of the given decomposition and check that each 
one of these leads to a consistent asymptotic scheme as in [3]. 

The equilibria of (|3. 14113.13)) correspond to the 'equilibria at infinity' of the flow of (|3.13p on the Poincare 
sphere <S2_, and they are precisely the points where 

IQsi - YP s i = (3.18) 
X 2 + Y 2 = 1, (3.19) 

with (P51 , Qgi)(X, Y) being the circle vector field introduced above. Therefore the study of the flow of the 
original system (|3.3p near its finite-time singularities is equivalent to that of the behaviour of the orbits of 
the asymptotic system (|3. 14113.15]) near its equilibria given by (|3.18H3.19|) . Because of the circle condition 
(|3.19p . these equilibria are of the form (X,Y, 0) and come in as antipodal pairs distributed along the circle 
of infinity, hence it is only necessary to study the flow along the equilibria having X > 0, and secondly those 
having Y > (in particular, they cannot both be zero). The flow will be topologically equivalent (with 
direction reversed) at the antipodal points if M, defined in Q3.7p . is odd (even) as in the homogeneous case, 

cf. 0E|. 

To study the flow near such equilibria, we follow a 'fan-out' technique, described in Rcfs. [7] for the 
all terms dominant case. Suppose we have an equilibrium, solution of ()3.18ti3.l5|i . with Y > 0. We define 
the fan- out map 

e = c = f (3- 2 °) 

which spreads any small neighborhood of that equilibrium on the circle back out onto the tangent plane 
to the Poincare sphere at Y = +1. Then the geometric dynamics defined by the asymptotic system (|3.14t - 
13.15)1 near this equilibrium is transformed into one described by the fan-out system near its corresponding 
equilibrium in terms of the new variables £, £: 
Fan-out system for Y > 0: 

-Mr, f £ M t-/-Mr\ 1 



C = =C M+1 Qsi(il)- (3-22) 



This is valid in a neighborhood of the point (0,+l,0). Similarly, there is a completely analogous fan-out 
reduction for any equilibrium having X > 0. The final result is: 
Fan-out system for X > 0: 

,Jfn ( 1 'A .aM n f 1 V 



* - ^(viJ-^'MviJ (3 ' 23) 

A = =-A»+'Q sl Q,|V (3.24) 

where the fan-out map in this case is rj = Y/X, A = Z/X. Note that in both fan-out systems we use the 
asymptotic decompositions of the field. 

We conclude that the dynamics of the original system near its singularities is determined by the nature 
of the equilibria of the corresponding fan-out systems given above, and these can be studied using standard 
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dynamical systems methods. In particular, if these turn out to be hyperbolic, then the situation will result 
in a relatively easy problem, whereas when we have non-hyperbolic equilibria, the dynamics in phase space 
'at infinity' will be more involved and the resulting global phase portraits will look more complicated. 

4 Discussion 

In this paper we have described a combination of the techniques of asymptotic splittings and central ex- 
tensions and provided an asymptotic from of the Poincare compactification suitable for the analysis of 
finite-time singularities. This applies to weight-homogeneous vector fields and their decompositions near a 
time when (some component of) the solution blows up. This description will be useful in treating problems 
in mathematical cosmology and other fields where such situations naturally arise. 

Although wc have described the asymptotics of the Poincare central extension method in dimensions 1 and 
2, there is nothing to prevent these results from being valid in any finite dimension. The 3-dimensional case 
of a homogeneous vector field with a singularity as t —¥ oo is described in Ref. [5] , and it is straightforward 
to formulate the results of the present paper in this context. In practise, however, we do not expect to deal 
with problems in dimensions higher than four, except in rare circumstances. For an abstract formulation in 
general dimensions on manifolds, we refer to [12] . 

The study of infinity in the context of the present paper opens the way for related structural stability stud- 
ies and applications to mathematical cosmology where such issues have especially important physical and ge- 
ometric interpretations, for instance a characterization of instabilities through the possible non-hyperbolicity 
of equilibria at infinity. Although, Peixoto's theorem for structurally stable vector fields concerns compact 
manifolds, there is a suitable extension of structural stability to the non-compact case, (for a modern intro- 
duction to these results, cf. Ref. [S]. 

Such results are of great importance in the present context, for although having the field blowing up 
only as t — > oo leads possibly to compact phase spaces and so stability may be natural for such vector fields 
when they are further compactificd on the sphere, the presence of finite-time singularities and the related 
possibility for orbits to escape to infinity in a finite time governed only by the dominant part of the vector 
field, opens real possibilities for the influence of non-compact phase spaces. This is eventually related to 
the possible validity of the index theorem under the precise conditions of the problem analyzed here, and 
it is perhaps interesting that known results in this direction, cf. [13] , are not directly relevant to the issues 
considered in this paper. 
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